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Summary 

Capacitance fuel gauges have served as the ba- 
sis for fuel quantity indicating systems in aircraft for 
several decades. However, there have been persistent 
reports by the airlines that these gauges often give 
faulty indications due to microbial growth and other 
contaminants in the fuel tanks. This report describes 
the results of a feasibility study of using gamma ray 
attenuation as the basis for measuring fuel quantity 
in the tanks. Studies with a weak Am 241 59.5-keV 
radiation source indicate that it is possible to contin- 
uously monitor the fuel quantity in the tanks to an 
accuracy of better than 1 percent. These measure- 
ments also indicate that there are easily measurable 
differences in the physical properties and resultant 
attenuation characteristics of JP-4, JP-5, and Jet A 
fuels. The experimental results, along with a sug- 
gested source-detector geometrical configuration, are 
described. 

Introduction 

Capacitance fuel gauges have served as the ba- 
sis for fuel quantity indicating systems in aircraft for 
several decades. These gauges, in the form of concen- 
tric cylinders, are mounted vertically at several loca- 
tions inside the fuel tanks (ref. 1). The summations 
of their indications give the total tank fuel content 
at any time. However, there have been persistent 
reports (ref. 2) by the airlines that the capacitance 
gauges often give faulty indications of tank fuel con- 
tents. The problem has been attributed to microbial 
growth and/or contaminants in the fuel tanks. The 
microbes can occur in storage tanks, delivery lines, 
pump trucks, and consequently, in the aircraft fuel 
tanks. The microbes attack the capacitance cylin- 
der coatings and thus expose the cylinder surfaces 
(electrodes) for subsequent corrosion and electrical 
noise in the capacitance bridge circuit. They also 
corrode the output signal leads. It is thus highly 
desirable that a fuel quantity indicating system in- 
sensitive to fuel contamination be developed. Such a 
system should be highly accurate (better than 1 per- 
cent), safe to use and operate, and inexpensive. 

An investigation of the feasibility of using gamma 
ray attenuation as the basis for measuring the fuel 
quantity in aircraft tanks has been conducted. The 
results of these studies are described in the following 
sections. 

Principle of Operation of a Nuclear Gauge 

The operation of a nuclear gauge is based on the 
attenuation of gamma rays passing through matter. 
As a result of interaction of gamma rays with the 
atoms in the test medium, the number of unaffected 


primary photons arriving at the detector is a func- 
tion of the path length in the test medium. For a 
uniform medium, it is given by the following expres- 
sion (ref. 3): 

I x *= l 0 e~i* ( 1 ) 

where 

l x number of unaffected primary 

photons transmitted through test 
medium 

I 0 number of photons incident on test 

medium 

/i linear attenuation coefficient for 

incident photons in test medium 

x path length in test medium 

Clearly, such a gauge will be more sensitive if the 
attenuation coefficient (//) is large for the incident 
photons. This dictates the choice of low-energy 
(less than 100 keV) photon sources. Two plausible 
candidate sources that meet the necessary criteria of 
low photon energy, long source half-life, and a well- 
resolved photon spectrum are Am 241 (458 years) and 
Cd 109 (453 days). 

The decay schemes (ref. 4) for these two sources 
are shown in figures 1 and 2, respectively. It is noted 
that the 59.5-keV radiation from the Am 241 source 
results from a super-allowed electric dipole (El) tran- 
sition in Np 237 , whereas the 87.7-keV radiation from 
the Cd 109 source arises from a weakly allowed electric 
octopole (E3) transition in Ag 109 . The latter tran- 
sition is strongly internally converted and produces 
a large, lower energy Ag K X-ray flux. For example, 
a 10-mCi Am 241 disc source emits 7.4 X 10 6 pho- 
tons (59.5 keV) per second per steradian, whereas 
a 10-mCi Cd 1( ^ 9 source emits 2.6 X 10 7 photons 
(22.6 keV) per second per steradian (ref. 5). The 
relative intensities of gamma rays and characteristic 
X-Rays emitted from these sources are summarized 
in table I (ref. 4). Thus, even though the choice of 
a Cd 109 source will necessitate changing the source 
every 3 years or so, Cd 109 still appears to be a viable 
candidate source by virtue of its large lower energy 
photon yield. 

Experimental Procedures for Measuring 
Attenuation Coefficients 

Since the exact compositions of aviation fuels are 
seldom known (refs. 6 and 7), it was not possible 
to calculate their attenuation coefficients for Am 241 
and Cd 109 gamma rays. It was therefore decided 
to determine the attenuation coefficients of selected 
types of fuels experimentally. 



Attenuation coefficients of several samples of com- 
mercial aviation fuels were measured in the narrow 
beam geometry illustrated in figure 3. The fuel cells 
were made of glass and were fabricated in the form of 
3-in. (7.62-cm) diameter flat-ended cylinders of three 
different lengths for easy data reduction. The gamma 
rays were detected with a 2-in. (5.08-cm) diameter x 
2-in. (5.08-cm) thick Nal (Tl) crystal coupled to a 
high-gain photomultiplier. Figure 4 shows the geo- 
metrical details of the source, collimators, fuel cells, 
and detector assembly. 

Measurements were made with empty fuel cells 
and cells filled with the test fluids. To further test 
the sensitivity of the system, measurements were also 
made with distilled water in the fuel cells. Typical 
Am 241 and Cd 109 spectra are shown in figures 5 
and 6. 

For the Am 241 source, the single-channel analyzer 
(SC A) limits were adjusted to accept the strong 59.5- 
keV peak. For the Cd 109 source, the SCA limits 
were set to accept the weaker 87.7-keV total capture 
peak rather than the stronger, but unresolved, lower 
energy Ag K a (22.1 keV) and Ag (25.0 keV) 
peaks. 

The nominal radioactive source strengths readily 
available for this test were of the order of 10 yzCi 
(Am 241 ) and 100 //Ci (Cd 109 ). They provided good 
counting statistics for all test fluids over a period 
of 10 minutes. Measurements were made with and 
without the source in each case to subtract the counts 
due to cosmic rays and other background sources of 
radiation. 

The attenuation coefficients were measured for 
water, JP-4 fuel, JP-5 fuel, Jet A fuel, regular 
leaded automobile gasoline, and unleaded automo- 
bile gasoline. 

Data Reduction and Results 

Counts were recorded for 10 minutes for each 
source for the three fuel cells filled with the test flu- 
ids. The geometrical details of the configurations in- 
corporating test ceils G-2, G-3, and G-4 are shown in 
figure 4. Typical results are summarized in table II. 

As illustrated in figure 4, the photons have to pass 
through air, glass fuel cell ends, test fluid, and a 
0.079-cm-thick aluminum housing for the Nal (Tl) 
crystal before arriving at the detector surface, that 
is, 

J x = J 0 ^g~^air x air e ~ ^ glass ^ glass 

X £‘~^ftuid* c fluid f J, Al x A\ ^ ^2) 

The values of /i a i r and pai at 59.5 keV and 87.7 keV 
have been reported by a number of authors (refs. 8 


to 10). With these values, I Q can be easily calculated 
from equation (2) if the entire path length is made 
up of air and aluminum. If an empty glass fuel cell 
is introduced in the path of the beam, the drop 
in the counting rate provides a direct measure of 
//glass for the incident photons. If the fuel cells 
are filled with the test fluids, the changes in the 
counting rates will reflect the effects of attenuation 
characteristics of the test fluids. The experimental 
values of linear attenuation coefficients of the various 
test fluids are summarized in table III. These values 
are based on several independent sets of data of the 
type summarized in table II. 

Since the mass attenuation coefficients 1 of the me- 
dia are of more fundamental value than the linear 
attenuation coefficients (ref. 3), it was necessary to 
determine the densities of the test fluids. The densi- 
ties of all the fluids were measured with a standard 
50-ml pycnometer, and these values were used to cal- 
culate the mass attenuation coefficients of the test 
fluids. These results are also included in table III. It 
is interesting to note that there are easily measurable 
differences in the attenuation coefficients of various 
test fluids. 

Subsequent to the measurements of the respective 
attenuation coefficients of all the test fluids for the 
Am 241 and Cd 109 gamma rays, it was finally decided 
to test the sensitivity of attenuation of low-energy 
photons as the basis for a fuel gauging system aboard 
aircraft. The wing tank geometry for a Boeing 737 
airplane was selected for the computer model as 
representative of all aircraft with positive wing tip 
inclination with respect to the horizontal while on 
the ground. The computational procedures and the 
program developed for calculations for an Am 241 
gamma ray source and Jet A fuel test medium are 
described below. 

Computational Procedure 

The computer program WNGTNK is written in 
FORTRAN Version 5 language for the Control Data 
CYBER 170 series digital computer system with 
network operating system (NOS) 2.3. The program 
requires approximately 40 000 octal locations of core 
storage. A typical case requires less than 4 central 
processing unit (CPU) seconds on the CYBER 173. 

The wing tank (Boeing 737) modeled by the pro- 
gram is illustrated in figure 7. For purposes of model- 
ing, each of the 14 compartments in this figure is ap- 
proximated by a rectangular box. Any similar wing 

1 The mass attenuation coefficients are independent of the 
actual density and physical state (gas, liquid, or solid) of the 
absorber. 
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tank can be modeled by this technique by simply ad- 
justing the number of compartments and the dimen- 
sions of each rectangular box. Figure 8 illustrates the 
tank model as viewed from the front of the aircraft 
with the fuselage (not shown) to the left. The pro- 
gram provides the user with the capability of speci- 
fying the height of the bottom of each compartment, 
B { , to simulate the bending of the wing associated 
with flight conditions. The solid dot (•) in each com- 
partment depicts the source location, SL, and the 
detectors are assumed to be fixed to the bottom of 
each compartment. Table IV summarizes the specific 
data used in the modeling of the Boeing 737, where 
Wj, Hi, and Pi are the compartment widths, heights, 
and depths, respectively. 

Once the tank geometry has been defined, the 
program steps through fixed percentages of tank fuel 
capacity. For each amount of fuel, the fuel level is 
computed with the assumption of a level fuel surface. 
With the fuel level known, the path length between 
each source-detector pair occupied by fuel or air is 
determined. From these path lengths, the number of 
counts is determined. The baffles between compart- 
ments contained in the wing structure are assumed to 
absorb radiation, so there is no interference between 
adjacent compartments. 

Program input consists of 14 numbers, separated 
by commas, representing the height, JB^, of each 
compartment bottom in inches. Program output 
includes both tabular and .graphic results. 

Typical results corresponding to the configuration 
of figure 8 are included as table V and are illustrated 
in figure 9. These data were acquired with a source 
strength of about 30 /iCi at each station in a count- 
ing interval of 1 second. Obviously, this system has a 
fast response time (approximately 1 second) and high 
resolution (approximately 1 percent). In this figure, 
each line depicts the relationship between counts and 
fuel expended for a specific compartment, with the 
lines toward the right nearer the wing tip and the 
lines toward the left nearer the fuselage. In partic- 
ular, note that when the tank is full, the counting 
rates are the same in each compartment, since the 
path lengths through fuel are all equal. As fuel is ex- 
pended, the counting rates change first in those com- 
partments near the wing tip. After approximately 
35 percent of the fuel has been expended, the com- 
partment nearest the tip is empty and shows no fur- 
ther change in counting rate. Also note that the 
source in compartment 1 is completely immersed in 
the fuel until approximately 85 percent of the fuel is 
expended and begins to show a change in counting 
rate as the fuel is reduced below this level. Figure 9 
also shows that significant changes in counts can be 


observed in one or more compartments as the fuel 
level varies, regardless of the tank contents, 

A listing of the computer program used in this 
analysis is included as an appendix. 

Discussion 

For the sake of specificity, we will confine our 
discussions to the results for an Am 241 (59.5-keV) 
gamma source. Similar results are expected for a 
Cd 109 (87.7-keV) gamma source. 

As seen from the data in table V, the counting 
rate is constant at all stations when the tank is full. 
A 1-percent reduction in the fuel content in the tank 
causes a large increase (about 56.9 percent) in the 
counting rate at the wing tip detector (station 14). 
A further reduction of 1 percent in the fuel causes 
an additional increase (about 26.6 percent) in the 
counting rate at the wing tip detector. It also results 
in a counting rate increase of about 16.2 percent at 
station 13. These counting rate changes are easily 
measurable. The same trend continues as more fuel is 
consumed. For example, when 10 percent of the fuel 
has been consumed, the total cumulative counting 
rate increases at stations 14, 13, and 12 are 238.6, 
114.6, and 35.9 percent, respectively. At the other 
end of the spectrum when the tank is nearly empty, 
the counting rates in the outer station detectors have 
stabilized, but the counting rates at the stations 
near the fuselage are changing fast. For example, 
when the tank is only 5 percent full, the counting 
rates at stations 1, 2, and 3 are 330.1, 521.0, and 
822.6 percent higher than the counting rate for the 
full tank. A further reduction of 1 percent in the 
fuel causes the counting rates to increase to 392.3, 
619.1, and 977.1 percent of the values for the full 
tank, respectively. • 

From these data it is apparent that the fuel 
quantity gauging system detailed in this report is 
capable of detecting changes as low as 1 percent in 
the fuel contents at the two extreme limits, that is, 
when the tank is almost full and when it is almost 
empty. A careful examination of table V illustrates 
that a similar degree of sensitivisty exists for all levels 
of tank fuel contents. 

From the foregoing discussion, it is apparent that 
a continuous monitoring of counting rates at all the 
detector stations should enable continuous tracking 
of airplane fuel tank contents with a high degree of 
sensitivity. 

Concluding Remarks 

It has been demonstrated that a suitably designed 
nuclear gauge should enable a continuous monitoring 
of the tank fuel contents to an accuracy of better 
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than 1 percent. Such accurate information — both 
at the point of flight origination when the tanks are 
presumably full and at the final destination when the 
tanks are almost empty — should prove very useful to 
the airlines. It should provide reliable information 
about the payload capacity at the beginning of the 
flight and safety margin near the end of the flight. 
The nuclear gauge is not expected to be susceptible 
to the fouling and corrosion problems experienced by 
the conventional capacitance gauges, since both the 
source and the radiation detector are sealed. Any 
algae or microbial growth on the source and detector 
windows can be easily removed during scheduled 
periodic maintenance checks of the gauging system. 

An added advantage of the nuclear gauge is its in- 
herent capability to detect water buildup in the tank. 
Since water is expected to gravitate toward the fuse- 
lage, any reduction in the counting rates at stations 1 
through 5 when the tank is at least half full can be 
used to infer the quantity of water in the tank. It 
is also a self-calibrating system with a high degree 
of cross-checking capability. This capability renders 
the nuclear gauging system independent of any back- 
ground count rate changes with altitude. (In any 
case, changes in background count rate at altitudes 
less than 10 miles are expected to be minimal in the 
SCA window centered at 59.5 keV.) 

It should perhaps be noted that despite the large 
low-energy photon flux obtainable with a Cd 109 
source, an Am 241 source would be more economical, 
since it would require no source replacement because 
of its long half-life. It would also be comparatively 
safer to handle and/or shield because of its lower en- 
ergy. As a matter of fact, Am 241 -based densitome- 
ters are currently in use aboard some aircraft. The 
licensing requirements for an Am 241 -based fuel quan- 
tity measurement system would be no different from 
what they are for those aircraft. By an appropriate 
choice of the Am 241 source strength, the response 
time of the nuclear gauge can be safely arranged to 
be less than 1 second. 


The effects of temperature on the fuel vol- 
ume can be easily taken care of by simultaneous — 
but independent— measurements of temperature and 
density. These measurements will also enable real- 
time computation of fuel mass (as opposed to fuel 
volume) at any time in flight or on the ground. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 6, 1986 
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TABLE I. RELATIVE INTENSITIES OF CHARACTERISTIC 
X-RAYS AND GAMMA RAYS EMITTED FROM Am 241 
AND Cd 109 RADIOACTIVE SOURCES 


Am' 241 source 

Cd lu ^ source 

Photon 


Photon 


energy 

Relative 

energy 

Relative 

(keV) 

intensity 

(keV) 

intensity 

11.89 (Np Li) 

2.2 

22.1 (Ag Ka) 

25.5 

13.90 (Np L a ) 

37.5 

25.0 (Ag Kp) 

5.0 

17.80 (Np Lp) 

51.2 


20.80 (Np L 7 ) 

13.8 

87.7 : | 

1.0 

26.35 

7.0 

| 


59.50 

100.0 




TABLE II. COUNTS PER 10-MINUTE INTERVAL FOR VARIOUS TEST MEDIA 
WITH Am 241 AND Cd 109 SOURCES 3 - 


Test 

medium 

Am 241 source 

Cd 109 source 

Cell G-2 
1*2 = 4.982 cm 

Cell G-3 
1.2 = 7.522 cm 

Cell G-4 
= 10.062 cm 

Cell G-2 
^2 — 4.982 cm 

Cell G-3 
L 2 = 7.522 cm 

Cell G-4 
£2 — 10.062 cm 

Air (no cell) 

75446 

49988 

35710 

43650 

30476 

23184 

Air (empty cell) 

54489 

36705 

25491 

35404 

24897 

18906 

JP-4 fuel 

27880 

13899 

7480 

21080 

12950 

9192 

JP-5 fuel 

26859 

13204 

7018 

20765 

12569 

9018 

Jet A fuel 

26816 

12918 

7229 

20322 

12431 

8802 

Leaded gasoline 

28732 

14462 

7981 

21110 

13105 

9233 

Unleaded gasoline 

28966 

14353 

8094 

20948 

13107 

9202 

Water 

22095 

9832 

5267 

18455 

10937 

8001 

Background 

1817 

1850 , 

1876 

5319 

5391 

5270 


a See figure 4 for geometrical details of fuel cell and associated shields/ collimators. 
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TABLE III. SUMMARY OF ATTENUATION COEFFICIENTS FOR VARIOUS TEST FLUIDS 


Test fluid 

Test fluid 
density, p, 
g/cm 3 

Am 241 (59.5 keV) source 

Cd i(Jy (87.7 

ieV) source 

p, cm 1 

Hm cm 2 /g 

p, cm -1 

fi m , cm 2 /g 

JP-4 fuel 

0.7546 

0.143 ± 0.003 

0.190 ± 0.004 

0.127 ± 0.002 

0.169 ± 0.003 

JP-5 fuel 

0.8097 

0.150 ± 0.002 

0.185 ± 0.003 

0.134 ± 0.004 

0.165 ± 0.005 

Jet A fuel 

0.8107 

0.150 ± 0.002 

0.185 ± 0.003 

0.137 ± 0.002 

0.168 ± 0.003 

Leaded gasoline 

0.7300 

0.135 ± 0.001 

0.185 ± 0.002 

0.126 ± 0.003 

0.172 ± 0.004 

Unleaded gasoline 

0.7443 

0.135 ± 0.002 

0.182 ± 0.003 

0.125 ± 0.002 

0.167 ± 0.003 

Water 

0.9974 

0.194 ± 0.002 

0.194 ± 0.002 

0.165 ± 0.002 

0.166 ± 0.002 


TABLE IV. DATA USED FOR BOEING 737 WING TANK MODEL 

Source type — Am 241 (59.5 keV); source strength — 10 6 counts per second; 
source enclosure — O.Ol-in.-thick aluminum 


Compartment 

w, 

in. 

H, 

in. 

D, , 

in. 

B, 

in. 

SL, 

in. 

i 

24.0 

26.8 

82.0 

0.0 

8.2 

2 

24.0 

24.6 

78.0 

1.5 

8.2 

3 

24.0 

22.6 

73.0 

3.0 

8.2 

4 

24.0 

20.6 

68.0 

4.5 

8.2 

5 

24.0 

18.6 

63.0 

6.0 

8.2 

6 

24.0 

16.6 

58.0 

7.5 

8.2 

7 

24.0 

14.6 

53.0 

9.0 

8.2 

8 

24.0 

13.6 

49.8 

10.5 

8.2 

9 

24.0 

12.7 

46.5 

12.0 

8.2 

10 

24.0 

11.8 

43.2 

13.5 

8.2 

11 

24.0 

10.9 

39.9 

15.0 

8.2 

12 

24.0 

10.0 

36.6 

16.5 

8.2 

13 

24.0 

9.1 

33.3 

18.0 

8.2 

14 

24.0 

8.2 

30.0 

19.5 

8.2 
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Appendix 

Listing of Computer Programs 
Program WNGTNK 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
56 


PROGRAM WNGTNK ( OUTPUT, INPUT, TAPE6-0UTPUT, TAPES- INPUT) 
CnMM0N/GE0MTY/W(14),H<14),D<l4),B(14)»SL(14)>VTaT,NC 
COMMON/ANS /COUNTS (14,101) 

DIMENSION FRAC(lOl) 

INTEGER COUNTS 
DATA ERROR/1. E-5/ 

1 F0RMAT(1H1,45X, ’SUMMARY OF COUNTS*//// 

1 5X,'* LIQUID * , 45X, 'COMPARTMENT' / 

2 18X, * 1 » 2 3 4 ',* 5 », 

3 * 6 7 • 8 *,« 9 «,* 10 *, 

4 * 11 12 »»• 13 14 '//) 

2 F0RMAT(1H0,5X,F7.2,2X,14I7) 

CALL PSEUDO 

CALL CALPLT{1.,1«,— 3) 

CALL IN IT 
NCM - NC - 1 
NFILL - 101 
READ 1 5,*, END- 10 ) B 
10 IF ( EOF( 5) .NF.O) GO TO 120 
ISGN » -1 

V - -0.01*VT1T 
DO 100 K-l, NFILL 
ISGN » -ISGN 
DVOL « 0.01*VT0T 

V - V + DVOL 
FRAC(K) - 100.+V/VT0T 
7-0. 

DO 20 I -1, NCM 
7 - 8(1+1) 

VI - 0. 

DO 20 J-1,I 
TOP = Z 
BOTTOM - B ( 4 ) 

IFCTOP.GT. (B(J) + H ( 4 ) ) ) TOP « B(J) + H(J) 

VI - VI + D(J )*W ( J )* (TOP - BOTTOM) ■ 

IF(Vl.GT.V) GO TO 30 
20 COMTINUF 
I - NC 
30 CONTINUE 
J - I - 1 
VPOT - 0. 

IF(J.fcO.O) GO TO 50 
DO 40 I-1,J 
Z - B ( J +1 ) 

TOP « Z 
BOTTOM « B (I) 

I F ( TOP , GT . (B( I) + H ( I ) ) ) TOP « B(I) + HII) 

V BOT « VBOT + D ( I ) *W ( I ) *( TOP - BOTTOM) 

40 CONTINUE 
50 CONTINUE 
4-4 + 1 
70LD - 8 ( J ) 

07 - Bt 4 ) + 0.7 5*H ( 4 ) 

60 CONTINUE 
VI - 0. 

ZM AX - 0. 

DO 70 1-1,4 

IF( (B (I ) + H(I) J.GT.ZMAX) ZMAX - B(I) + H(I) 
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59 


TOO ■ ZOLD 

60 


BOTTOM * B ( J ) 

61 


IF ( TOP.GT .(8(1) + H ( I ) ) ) TOP * B(I 

62 


IF (BOTTOM. GT. (B (I) + H(I ) ) > BOTTOM 

63 


VI • VI + D(I)*W(I)*(TOP - BOTTOM) 

64 


70 CONTINUE 

65 


VOLUME - VBOT + VI 

66 


TST V * V 

67 


IFIV.EQ.O.) TSTV - 1. 

68 


IF(ABS( VOLUME - V ) /TS TV .LT . ERROR ) ' 

69 


IFCVOLUME.LT. V) GO TO 80 

70 


TOLD * ZOLD - 07 

71 

■- 

07 * 0.5*07 

72 


GO TO 60 

73 


80 ZOLD = ZOLD + 02 

74 


GO TO 60 

75 


QO CONTINUE 

76 


IF(ZOLD.GT.ZMAX) ZOLD « ZMAX 

77 


CALL TABL ( ZOLD» K) 

78 


IF CK.EQ.l ) CALL PICT(ZOLO) 

79 


100 CONTINUE 

80 


WPITE(6>1) 

81 


DO 110 I«1,NFILL 

82 


WR ITE (6> 2 ) FR AC (I), (COUNTS (K#I),K- 

83 


110 CONTINUE 

84 


CALL PLTCNT 

85 


120 CONTINUE 

86 


CALL CALPLT(0,>0.#R99) 

87 


STOP 

88 


END 


+ H( I ) 

* B ( I ) + 


H C I ) 


GO TO 90 


1 *NC ) 


Subroutine INIT 

1 

2 

3 

4 

5 

6 
7 
fi 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


SUBROUTINE INIT 

C0MM0N/GE0MTY/W(14)*H(14)>D(14)>B(14)>SL( 14)»VT0T»NC 


NC 
D ( 1 ) 

D ( 2 ) 
0(3) 
0(4) 

D (5 ) 

D ( 6 ) 
0(7) 
0 ( 8 ) 

D ( 9 ) 
0 ( 10 ) 
0 ( 11 ) 
D ( 12 ) 
D (13 ) 
0(14) 
H ( 1 ) 

H ( 2 ) 

H (3 ) 

H ( 4 ) 

H ( 5 ) 

H (6 ) 

H (7) 

H ( 8 ) 

H ( 9 ) 

H ( 10 ) 
H ( 11 ) 
H ( 12 ) 
H( 13 ) 
H ( 14 ) 
VTOT 
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82. 

78. 

73. 

68 . 

63. 

58. 

53. 

49.8 

46.5 

■ 43.2 

• 39.9 

■ 36.6 

* 33.3 
» 30. 

26.8 

24.6 

22.6 
20.6 
18.6 
16.6 

14.6 

13.6 

12.7 

« 11.8 

■ 10.9 

■ 10 . 

« 9.1 

« 8.2 

0, 



33 


DO 10 I *1» NC 

34 


W ( I ) « 24. 

35 


VTOT * VTOT 

36 

10 

CONTINUE 

37 


NC02 - NC/2 

38 


m 20 i * i»nc 

39 


SL ( I ) - 8.2 

40 

20 

CONTINUE 

41 


RETURN 

42 


END 


+ wm*HtI}*D(I) 


Subroutine PICT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
2P 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


SUBROUTINE PICT(Z) 

C0MM0N/GE0MTY/W(14)»H(14),D(14)#B(14),SL(14), VT0T>NC 
DIMENSION XA(35),YA(35),XB(35),YB(35) 

DATA EPS/1. €-4/ 

IND * 1 

X * 0. 

Y « 0. 

XI « 0. 

YI « 1 .OM 8 (NC ) + H(NC ) ) /34. 

CALL C A L P L T ( X I # Y I >3 ) 

XA(IND) - XI 
YA(IND) « YI 

IF(1.*Z/34..LT.YI) YA(IND) » 1.+Z/34. 
IF(YA(IND).LT.l»*B(NC)/34.) YA(IND) » l.*B(NC)/34. 

IND * IND + 1 

DD 10 I-l/NC 

X « X + W ( NC+l-I ) /34 • 

Y - 1 .OMB (NC+.l-I) + H(NC+1-I) )/34. 

CALL CALPLT(X,Y,2) 

XA(IND) » X - EPS 
YA(IND) * Y 

IF(1.*Z/34,.LT.Y) YA(IND) » l.M/34. 

IF ( YA ( IND) .LT.l »*B (NC+l-I ) /34. ) YA(IND) « 1 . *B ( NC+l-I) /34. 

IND “ IND + 1 

IF(I.NE.NC) 

1Y » 1 .0* ( B (NC — I ) + H( NC— I ) )/34. 

IF(I.EQ.NC) Y » 0. 

CALL CALPLT(X#Y,2) 

XA(IND) * X + EPS 
YA ( IND) » Y 

IF ( 1 «*Z/34 . .LT. Y) YA(INO) « l.*Z/34. 

IF(YAdND) .LT,l.*B(NC-I)/:34.) YA ( IND) - 1 .*B ( NC-I ) /34 . 

IND =* IND + 1 
10 CONTINUE 

CALL CAL PL T(XI»YI>3) 

X - XI 

Y * l,0*B(NC)/34. 

CALL CALPLT (X»Y>2) 

IND - 1 

XB(IND) * X 

YB(IND) - Y 

IND * IND + 1 

DO 20 I « 1 t NC 

X - X + W (NC+l-I ) /34, 

Y - 1.0*B(NC+1— I)/34. 

CALL CALPLT (X> Y*2) 

XB(IND) » X - EPS 
YB(IND) - Y 

IND « INO + 1 
IF(I.EQ.NC) GO TO 20 
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51 


Y * 1.0*B(NC-I )/34. 

52 


CALL CALPLT(X,Y,2> 

53 


XB(IND) * X + EPS 

54 


YP(IND) « Y 

55 


IND * IND ♦ 1 

56 

20 

CONTINUE 

57 


XA(29) * 0. 

58 


X B ( 2 ° ) * 0. 

59 


X A ( 30 ) * 1. 

60 


X B ( 30 ) - 1. 

61 


YA ( 29 ) « 0. 

62 


YB ( 29 ) « 0. 

63 


YA (30 ) * 1. 

64 


YB (30 ) - 1, 

65 


NP « 28 

66 

C 

CALL HAFTONK XA,YA,NP,XB,YB,NP,9) 

67 


CALL LINPLT(XA,YA,NP,1,0,0,0,0) 

68 


XI * -0.5*W(14)/34» 

69 


DO 30 1*1, NC 

70 


XI * XI + W(NC-I+1> /34. 

71 


Y I - 1 »*(B (NC-I +1) + SL(NC-I+l))/34. 

72 


CALL PNTPLt (X I, YI, 22,1 ) 

73 

30 

CONTINUE 

74 


CALL NFRAME 

75 


CALL CALPLT{l.,l.,-3) 

76 


PFTURN 

77 


FND 


Subroutine TABL 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


SUBPOUT IN P TABL ( Z,K) 

C0MM0N/GE0MTY/W(14),H(14),D(14),B(14),SL(14),VT0T,NC 
COMMON/ AN S /COUNTS (14# 101 ) 

PEAL LIQSIG,LIQDEN 
INTEGER COUNTS 

DATA AIPSIG,AIRDEN,LIQSIG,LIQDEN,ALSIG, ALDEN/ 

1 0.178,0.001293,0.150, .7999,0.248,2.7/ 

DATA STRNGT/l.E+6/ 

PI * AC-OSC-1. ) 

DO 10 1*1, NC 
HEIGHT » H ( I ) 

TOP « Z 

IF ( TOP. GT » (H( I ) +B ( I ) ) ) TOP * H(I) ♦ B ( I ) 

IF (TOP. LT. B C I ) ) TOP * B(I) 

SLDANG * 10.*2.54/(4.*PI*(2.54*SL(in**2) 

DLIQ « TOP - BCI) 

DAIR « SL ( I ) - DLIQ 
IF (DLIQ.GT »SL Cl)) DLIQ * SL(I) 

IF ( DLIQ .EQ.SL ( I ) ) DAIR « 0. 

XC * SLDANG*S TRNGT 

XC • XC*£XP(-2.54*0.01*ALSIG*AL0EN) 

XC * XC*EXP(-2. 54*DLIO*LIOSIG*LIQDEN) 

XC * XC*EXP(-2.54*DAIR*AIPSIG*AIRD£N) 

COUNTS ( I ,K ) * XC 
10 CONTINUE 
RETURN 
END 



Subroutine PLTCNT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

1* 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


SUBROUTINE PLTCNT 

C0MM0N/GF0MTY/W(14),H<14),D(14)>B<14y,SL(14),VTnT>NC 
COMMON/ ANS /COUNTS (14,101) 

DIMENSION FRAC(lOl) 

INTEGER COUNTS 
DIMENSION X (105 ) * Y (105 ) 

DO 30 J-1,NC 

X ( 102 ) - 0. 

X ( 103) « 10. 

Y ( 102 ) * 0. 

Y (103 ) - 600. 

IFCJ.NE.l) GO TO 10 

CALL AXES ( 0«»0* *0.; 10. » X( 102) > X ( 103) # 1*> 0»# 

1 8H* LIQUIDS. 2*-8) 

CALL AXES (0.»0.»90./8.»Y{102),Y(103),1.#0.» 

1 6HC0UNTS j 0.2> 6 ) 

10 CONTINUE 

DO 20 1-1,101 
Xtl) * FLOAT 1 1-1 ) 

Y ( I ) • COUNTS ( J > I ) 

20 CONTINUE 

CALL LINPLT(X,Y, 101,1, 0,0, 0,0) 

30 CONTINUE 
RETURN 
END 
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Figure 3. Schematic diagram of experimental system used for measuring attenuation coefficients of Am 241 and 
Cd 109 gamma rays. 



Lj * 0.340cm (glass) 

1-2 * 4. 982 cm, 7. 522 cm, 10. 062 cm 
l _3 ■ 2. 540 cm (air) 

L 4 * 0. 079 cm (aluminum) 

S * Radioactive source 


G-2 * Glass cell with L 2 * 4. 982 cm 
G-3 * Glass cell with l_ 2 = 7. 322 cm 
G-4 * Glass cell with * 10. 062 cm 


Figure 4. Geometrical details of fuel cell and associated shields/collimators. 
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o 20 40 60 80 100 


Channel number 

(a) Empty cell. 



0 20 40 60 80 100 


Channel number 

(b) Cell filled with Jet A fuel. 

Figure 5. Typical spectra of Am 241 (59.5 keV) radiation source through measurement cell. 
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(b) Cell filled with Jet A fuel. 

Figure 6. Typical spectra of Cd 109 (87.7 keV) radiation source through measurement cell. 
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Fuselage 



Figure 7. Wing compartment diagram for Boeing 737 airplane. 
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Fuel content, percent , 

Figure 9. Counting rate versus fuel content in wing tank at various source-detector stations. 
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